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We analyze the tunnel coupling between an impurity state located in a (5-layer and the 2D de- 
localized states in the quantum well (QW) located at a few nanometers from the S - layer. The 
problem is formulated in terms of Anderson-Fano model as configuration interaction between the 
carrier bound state at the impurity and the continuum of delocalized states in the QW. An effect of 
this interaction on the interband optical transitions in the QW is analyzed. The results are discussed 
regarding the series of experiments on the GaAs structures with a iS-Mn layer. 

PACS numbers: 75.75.-c, 78.55. Cr, 78.67.De 

I. INTRODUCTION 

The problem of so-called configuration interaction of a single bound state with a continuum of states goes back to the 
famous paper by U. Fano^ rated as one of the most relevant works of 20th century;^ The suggested theoretical approach 
often regarded as Fano-Anderson model or configuration interaction succeeded in explaining puzzling asymmetric 
resonances observed in various experiments in atomic and nuclear physics, condensed matter physics and optics.—. 
The co-existence of the discrete energy level and the continuum states within the same energy range is also quite 
common in low-dimensional semiconductor structures 3"— Of particular interest nowadays are the structures having a 
quantum well (QW) and a ferromagnetic or paramagnetic layer located in the vicinity of the QW, but not penetrating 
into the QW region. In such structures high mobility of the carriers along the QW is combined with the magnetic 
properties provided by the magnetic layer. A number of recent experiments show that the Mn (5-layer gives rise to 
circular polarization of the photoluminescence (PL) from the QW in an external magnetic field applied perpendicular 
to the QW planci^ii It was questioned whether the spin polarization of the carries in the QW is due to the electrons 
tunneling to Mn site or the tunnel coupling of the holes at Mn with those in the QW. The latter mechanism seemed 
to lack the proper theoretical description. In this paper we try to fill this gap. We show that the simple scheme of 
the holes configuration interaction leads to the opposite sign of the circular polarization than that observed in the 
experiment. The model system considered in the present paper consists of a (5-layer of the impurities (donors or 
acceptors) and a QW having one level of size quantization for the electrons or holes respectively. The energy level 
of the impurity bound state lies within the range of the 2D states size quantization subband in the QW. We will 
be considering the case of rather deep impurity level in the sense that the impurity activation energy substantially 
exceeds the kinetic energy of the 2D carriers in the QW. The attracting potential of the impurity is assumed spherically 
symmetric and since it is a deep level we treat it with zero radius potential approximatiort^. At that we consider both 
the simple band structure and the one of the GaAs valence band type. 

II. TUNNELING BETWEEN IMPURITY AND QUANTUM WELL 

In this section we consider the configuration interaction between a single impurity bound state and the continuum of 
2D states in the QW. The potential barrier separating the impurity from the QW is assumed to be weakly transparent 
for the tunneling. Rigorous calculation of the eigenfunctions is rather hard to perform as it requires solving stationary 
Schrodinger equation in the complicated 3D potential. In order to circumvent the explicit solving of the Schrodingcr 
equation for tunneling problems the so-called tunneling or transfer Hamiltonian formalism is commonly used as 
originally proposed by Bardeen^. The total Hamiltonian is expressed as H = Hi + Hqw -\- Ht, where Hi is partial 
Hamiltonian having the bound state at the impurity as its eigen state. Hqw in the same way corresponds to the 
QW itself, its eigenfunctions ipx form non-degenerate continuum of states characterized by the quantum number (s) 
A. The term Ht accounts for the tunneling. In the secondary quantization representation the total Hamiltonian can 
be written as follows: 

H = eoa^a + J exc1^cxd\+ J (t\cla + t*xa^cx) dX, (1) 

where a^,a - the creation and annihilation operators for the bound state characterized by its energy Eq, and c^,ca 
- the creation and annihilation operators for a continuum state having energy ex- The energy here and below is 
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measured from the level of size quantization of the carriers in the QW so that e\ is simply their kinetic energy. The 
expression ([T} is rather general, in fact it can be regarded as introducing the coupling between two systems into the 
Hamiltonian in the most simple phenomenological way. From this viewpoint the coupling parameter t\ is still to be 
determined through exact solving of the eigenvalue problem for the whole system. Bardeen's approach suggests a 
simple recipe for calculation of the tunneling parameter for the case of weak tunneling through a potential barrier: 

tx= f i^lK^ - ^K^Ddv, (2) 

J a 

where integration is performed over region a to the one side of the barrier. Here K is the kinetic energy operator: 

K = -^A. (3) 
2m ^ ' 

The attraction potential of the impurity is considered spherically symmetric, so the whole system (impurity+QW) has 
the cylindrical symmetry with z axis directed normally to the QW plane and going through the impurity center. Thus 
for further calculations it will be most convenient to represent the QW states in cylindrical coordinates rather than 
as plane waves. In this case each state is characterized by the wavenumber k and the cylindrical harmonic number I: 



VH^v{z)^j:^Ji{kp)e^'' (4) 

where Ji{kp) is the Bessel function of order I, p and 6 are the polar coordinates in the QW plane, m-the in- 
plane effective mass, rj (z) is the envelope function of size quantization in z-direction. The wavefunction ^ has 
the normalization: 

{ipkil'Pk'i') = ~ (5) 

where e = h?k^ /2'm. The potential barrier separating the deep impurity level from the QW in the first approximation 
can be assumed having a rectangular shape. Inside the barrier the function r]{z) is (z-axis is directed towards the 
impurity, z = corresponds to the QW boundary): 

77(z)^^e-'^ (6) 
y a 



where 9 = y ; o, is the QW width, Eq is the binding energy of the bound state, at the same time Eq determines 
the height of the potential barrier. Let us firstly consider the simple band case valid for the bound electrons at donor 
impurity coupled to the QW conductance band. The spherical potential of the impurity results in the ground state of 
the carrier to be angular independent, therefore the efficient tunneling overlap occurs only with the zeroth cylindrical 
harmonic (pko = <p(e)- For the deep impurity level one can use zero radius potential approximation^ and express the 
s-type wavefunction as: 

i; = ^q—-. (7) 

The integration ^ over the space is reduced to the integration over the surface fls inside the barrier which is more 
convenient to take at the impurity site. This yields for the electrons tunneling between the donor state and the QW: 
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It is clearly seen that as long as the case k << q is considered, the tunneling parameter has very weak dependence 
on k. In order to apply the same approach to the holes tunneling in GaAs it has to be generalized for the case of the 
valence band complex structure. Let us consider In2;Gai_2.As QW having only one level of size quantization for the 
heavy holes and neglect the light holes being split off due to the size quantization. The basis of Bloch amplitudes to 
be used is formed of the states with certain projection of the total angular momentum J = 3/2 on z axis. It would be 
tempting to generalize ([2]) by treating K as the kinetic part of the Luttinger Hamiltonian [hk^fikyfikz are, as usual, 
the momentum operators along the appropriate axis): 
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The functions tpa, in © become now 4-coniponent vector functions (also the spin indices a and (3 are added 
here). The exphcit expression for the bound hole state functions Tpa and the 2D hole states i^^^ can be found in 
Ref J^. The important thing about those is while the decay length in z-direction of the 2D wavefunctions (p\p is 
controlled by the heavy hole mass ruhh ~ 0.5 mg (mg is the free electron mass), the decay length of radial part of the 
bound state wavcfunction ipa is characterized by both heavy hole mass rnhh and the light hole mass mni w 0.08 rrio^. 
Analogously to the simple band case the integration ^ over the whole space is reduced to the integration over the 
surface fls inside the barrier, at that, only z-projection of the kinetic energy operator is required. The expression for 
tunneling parameter simplifies into: 

tklc.p^''^ = {B- A) j^^ ds(^ipkW*^ija-iJa-^'fkW*y (11) 

where ipki is given by (|4]). 

Regrettably, the above given straightforward generalization of ([2]) fails to be fully correct. Indeed, the largest decay 
length of the bound state ipa is determined by the light hole mass while the decay length of the QW states is governed 
by the heavy hole. Due to this circumstance the result of the surface integration ([TT|) becomes dependent on the 
particular position of the integration surface inside the barrier. However, it can be shown that in the case of two 
masses the exponential dependence of the tunneling parameter on the barrier thickness is determined by the smallest 
mass, but the exact value of the tunneling parameter cannot be correctly obtained within the given approach. Now 
we define q = , /3 = mih/mhh- The explicit evaluation of the overlap integrals with account for /c << q 

shows that the tunneling configuration interaction to be accounted for is only between the zeroth cylindrical harmonic 
fko -| and the bound state as well as between ipi.Q ^3 and Both arc governed by the same tunneling 

parameter f^: 

where 1<X^2, C^l are weak dimensionlcss functions of k/q, m'f^f^ is the effective in-plane heavy hole mass. The 
tunneling parameter exponentially depends on the barrier thickness with the light hole mass entering the exponent 
index. The particular expressions for x and ( depend on the surface one chooses for the integration in ([2]). 

In both cases for t'j.,t'^ it is reasonable to assume that the tunneling parameter does not depend on k as weak 
tunneling implies k « q. Still, its rapidly decreasing behavior for k » q has to be kept in mind when it provides 
convergence for integration over k. In our estimations the shape of the potential barrier separating the QW was 
assumed rectangular. This is quite reasonable for the estimation at fc << g. However, the particular shape of 
the barrier becomes important when one is concerned with experimental dependence on the distance d between the 
impurity and the QW. 



III. EFFECT ON THE LUMINESCENCE SPECTRUM 



The transfer Hamiltonian ([T]) with known tunneling parameter t{e) allows one to construct the eigenfunctions of 
the whole system given those of the bound state V' and the QW states (p{e) : 

poo 

^{E)=vo{E)^+ v{E,e)^{e)de, (13) 
Jo 

E denotes the energy of the state 5*. Here ^p{s) are the wavefunctions with zeroth cylindrical harmonic, as was 
shown above the other harmonics are not affected by the tunneling configuration interaction. Plugging (|13p into the 
stationary Schrodinger equation: 



= E-^ 
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with H being the effective Hamiltonian ([T|) one gets the following system of equations: 

iyoiE)eo + t[e)y{E,e)de = Eua[E) , , . 

iy{E,e)e + t{e)i^o{E) = Eu{E,e). ^ ' 

In the present work we consider the case of the bound level energy lying within the range of the continuum: Eq >> . 
For this case the solution is obtained as shown ini: 

nH^{E) + {E-e^f' 
V {E, e) ^ (E) [P^^^ +Z{E)t {E)5 {E - e)^ , (15) 

where 

— TW) — ' 

f°° (e) 

F{E)= Pjj^.de, (16) 
J a [E-e) 

P stands for the principal value and Eq is the center of configuration resonance, which appears to be slightly shifted 
from Eq. 

SoiE)=Sa + FiE}. (17) 



Because of /c << g it is reasonable to put t =const everywhere, except for (|16p where decrease of t at i? — > oo is 
necessary for convergence of the integral. 

In order to analyze the influence of the configuration interaction on the luminescence spectra we have to calculate 
matrix element of operator M describing interband radiative transitions between the hybridized wavefunction 5'(i5) 
and wavefunction of 2D the carrier in the other band of the QW which we denote by £,k'i', here k' is the magnitude 
of the wavevector, I' is the number of cylindrical harmonic analogously to (|4]). If, for instance, one considers the 
acceptor- type impurity then ^{E) is the hybridized wavefunction of the 2D holes and ^k'l' is the wavefunction of the 
2D electrons in the QW. We assume that (a) there are no radiative transitions between the bound state wavefunction 
■0 and the 2D carrier wavefunction ^k'l' in the other band thus the matrix element for transitions from the bound 
state: 



k'l' 



M = 0, (18) 



(b) the interband radiative transitions between the free 2D states in the QW are direct. According to (|4]) the 
wavefunctions (p{e) and S,{s') corresponding to the zeroth harmonic in the cylindrical basis are: 



ae')=C{z)\I^Jo{k'p), (19) 



where 



_ , _ V2m'e' 

rv — ; . rC — - ; 



77(2), C(z) - the appropriate size quantization functions in z-direction, m, m' are the in-plane masses of the electrons 
and holes respectively if the donor-type impurity is considered and vice versa for the acceptor case. Without the 
tunnel coupling the matrix element for the direct optical transitions between the states ip{e) and ^(e') is given by: 



Mo (£,£') = (e(e') \m\ ^(e)) = {k ~ k') , (20) 

where Uk is the appropriate dipole matrix element for the Bloch amplitudes. According to the above mentioned 
considerations it is only this matrix element that is affected by the tunnel coupling, preserving the matrix elements 
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corresponding to the transitions between other than the zeroth cyhndrical harmonic. Denoting by M the modified 
matrix element for transitions between the states 'i>{E) and (e') with the further use of the Fano theorjii one obtains: 



M{E,e' f ^ Mo{E,e'f 



1 - 



We proceed further with the Fermi's Golden Rule for the transition probability: 



2 roo nOO 

W{nw)^— / \M{E,e')\^f' {e')f{E)5{E + e' + Eg-huj)dEde', 



(21) 



(22) 



where Eg - the QW bandgap, hjj - the energy of the radiated photon, /, /' - the energy distribution functions for 
the carriers in the hybridized and intact bands respectively. Substituting (^0)) and (dJ) into ([^ one should treat 
correctly the delta-function for the wavenumbers of the zeroth cylindrical harmonic. It can be shown that: 

where S is the area of the QW. Then we arrive at: 



^ = -^/(^-) 1 A , (23) 



where 

e' f„ -1 



f[E^) = f {a-'E^)fiE^), 

huj — En 



E^ = 



1 + a-i 



m + to' ' 



a = m' /m^ 



while for the all cylindrical harmonics altogether the unperturbed optical transition rate yields: 



(24) 



Wo {h.) = [j-^Sj . (25) 

The result (j23p obtained for a single impurity can be applied to an ensemble of impurities provided their interaction 
between each other is weak compared to the tunnel coupling with the QW. In this case the sample area S should be 
replaced with n~^, n being the sheet concentration of the impurities in the delta-layer. After normalization by the 
area of the QW from (|23|1 .([25 | we finally get the spectral density of the luminescence intensity: 

/ {hw, fo) - /o (fiw) 1 - a (e~o) V^^TT^T^ ^ ' (26) 

Y TT^f^ + [E^ - £o) / 

where 

h 

a(£o) = 



7r3/V2mfo(l + a^^)' 



27rM2m 
h (Sw) = ^3 ./ [E^ 
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IV. POLARIZATION OF THE SPECTRA 



It follows from (j26p that the bound state lying within the energy range of the continuum causes a dip in the 
luminescence spectra emitted from the QW. If then for any reason the bound state is split the luminescence spectra 
will show the appropriate number of the dips shifted by the splitting energy A. If one considers the splitting in the 
magnetic field applied along z each of the split sublevels is characterized by certain projection of spin and interacts 
with only one of the 2D carriers spin subbands characterized by the same projection of spin. Thus, for each of the two 
circular polarizations cr"'', a~ of the light emitted from the QW one would expect one dip, its spectral position being 
different for (t+ and in accordance with the splitting energy A. As an example let us consider the GaAs-based 
QW and 2D heavy holes interacting via the tunneling configuration interaction with the bound state at an acceptor. 
This case is shown schematically in Fig.l. The 2D holes with the projections of total angular momentum j = +3/2 
and 2 = —3/2 recombine emitting respectively right- (tr^) and left- (cr~) circularly polarized light. In section|ll]it was 
shown that the heavy holes with j = — 3/2(j = -1-3/2) interact basically with the bound states V'-S (V'+i)- An external 
magnetic field applied along z would cause Zeeman splitting of the bound state energy level Eq into = Eq -f- A/2 
and £q" ~ — The splitting 1\ = — may also originate from exchange interaction of the holes with 

spin-polarized acceptor ions. Let us refer to the case of Mn ions having positive g-factor {g « 3, see Ref.— ). The hole 
is coupled to Mn in antiferromagnetic way thus the level corresponds to j = —3/2 and to j = -1-3/2. As follows 
from ([T7)) . (P^ the difference in the positions of the resonances (dips) E'^ and iJj corresponding to the bound state 
sublevels and is given by: 



A = K+ - e: 



A In 1 




(27) 



Unless the positions of the resonances are too close to the band edge the last term in (P7)) can be neglected and 
A = A = e J — . With account for the energy distribution functions for the holes and electrons the shifted positions 
of the resonances lead to the difference in the luminescence intensity for the opposite circular polarizations. In the 
discussed example of the antiferromagnetic alignment of the hole the luminescence spectra /^(fiw, £o^), ^ (fii^,£o ) 
having the resonance positions at and e,^ correspond to the circular polarizations a~ and cr"*" respectively. As can 
be seen from (pS)) the difference in the resonance positions A = — Eq^ leads to the integral polarization of the 
spectra if the distribution function j(E) significantly varies in the vicinity of Eq. This is illustrated in Fig. 2. The 
functions and are shown by blue and red solid lines respectively. The integral polarization is naturally defined 
as: 



P = 



J I- {huj) d (huj) ~ J 1+ {hjj) d (hu) 

P{a+)-P{a^) Bg Bg 



P{a+) + P{a-) 



2 / Iq (hu;) d {hui) 



With use of ^ this yields: 



■0 



2 / / {E) dE 



f (E) dE 



(28) 



The slow varying functions f{E) and eQ{E) in the integral may be assumed as constants taken at Eq, , the tunneling 
parameter will be treated as a constant in the whole range of interest {E) = . 
Then treating the expression in brackets as delta-functions we obtain: 



^ne^ f (4) (4)"'^' - / (^o ) (^o ) 



p = - 



,-1/2 



23/2 



JfiE)dE 



(29) 



Note that for the considered example the polarization degree appears to be negative. The positive sign would have 
appeared if the ferromagnetic coupling between the acceptor ion and the hole had been assumed. 
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V. THE ELECTROSTATIC EFFECT 



Because of the tunneling involved in the polarization of the luminescence one might reasonably expect very strong 
dependence of the polarization degree on the distance d between the (5-layer and the QW (i.e. the thickness of 
the spacer). However, the purely exponential dependence of the polarization on the barrier thickness appears to 
be weakened due to the electrostatic effect shown in FiglT] and explained below. Let us for simplicity consider the 
electrons distribution function being nearly constant within the configuration resonances. The holes are considered to 
have Fermi distribution function characterized by the chemical potential fj, and the temperature T. In the absence of 
an external optical pumping the holes in the QW are in thermodynamic equilibrium with the acceptors in the (5-layer, 
therefore they have the same chemical potential. Under low pumping conditions the already large concentration of 
the holes in the QW is not strongly violated, so it is reasonably to assume that the quasi Fermi levels of the holes at 
the acceptors and in the QW coincide, it means that eq = fi. Strictly speaking, this is valid for a single bound level, if 
the level is split so that Eq — e'^ = A, one should probably assume Eq = ji. From (|29p we get the following simplified 
expression: 

^ ^ tanh— — (30) 



25/2^^^3/2 2kT 

As we will show below both t and fi contribute to the dependence of the integral polarization P on the spacer thickness 
d and the QW depth Uq. The holes in the QW provide an electrical charge density estimated as cr = eNfi, where e is 
the elementary charge, N is the 2D density of states. The positively charged plane of the QW and negatively charged 
(5-layer of partly ionized acceptors separated by a distance d produce an electric field 

F='-^, (31) 

E 

E being dielectric constant of the material. Due to the electric field F the valence band edge at position of Mn layer 
appears to be shifted from the valence band edge just outside of the QW hy F ■ d. Because the quasi Fermi level of 
the acceptors exceeds the local position of the valence band edge by the binding energy Eq, the equality of the quasi 
Fermi levels leads to a simple equation (see Fig[T]) : 

Uo = ^l + Eo + eFd, (32) 

where Uo is the QW depth and fi is the chemical potential of the holes in the QW. With (PT|) one gets : 



_ Ug-Eo {Uo-Eo)e 
^ ~ 1 + ^ 'inNed ■ ^ ' 



In order to estimate the dependence of the tunneling parameter t on the QW and spacer parameters we consider the 
WKB tunneling through trapezoid barrier as seen in FigH] With taking into account and ([55]) this leads to the 
following expression (we assume fi <ti Uq): 

f - exp {-Kd) , (34) 

where 

From (pO| . (p4)) . ((35)) . (|33)) follows the dependence of integral polarization on the spacer thickness: 

P ^ d^/^exp{-Kd), (36) 

Note that electrostatic effect results in the dependence of /i on d which leads to the dependence of P on c/ being 
not purely exponential but weakened by the pre-exponential factor d^/^. While the correction is pre-exponential it 
appears to be significant enough up to kc? ~ 2 — 3 which is typical for the experimental situation. 



VI. DISCUSSION 



In the proposed theory the polarization of light emitted from the QW originates from the splitting of the impurity 
bound state and therefore may exceed the polarization degree expected from an intrinsic g-factor of the 2D carriers 
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located in the QW. The sign of the polarization deserves special discussion. As was shown above, the tunnel coupling 
causes a dip in the luminescence spectra. This means that in the considered scheme the polarization of the luminescence 
from the QW is expected to be of the opposite sign than that due to the optical transitions between the bound state 
and the free carriers inside the barrier. In particular, the configuration interaction between the 2D heavy holes and 
Mn (5-layer considered in Sec lIVI leads to the negative sign of the polarization (a mistake made iniS has mislead to 
the positive sign). Such result contradicts the known experimental dataii^,— , where the polarization is shown to be 
positive. This might suggest that regarding these particular experiments the polarization is not due to the holes 
configuration interaction but rather due to polarization of the electrons as suggested ii>i^. The other possibility might 
be that the relevant bound state of the hole at Mn is more complex and docs not resemble the simple antiferromagnetic 
exchange coupling with Mn ion. 

Let us estimate the expected magnitude of the circular polarization degree due to the tunneling configuration 
interaction. We assume the deep impurity level Eq = 100 meV, the barrier thickness d = 5 nm, the QW width a = 10 
nm. Taking the effective mass as that of the electrons in GaAs m = 0.06 toq for the simple band case described by 
([HI) one gets for the tunneling parameter (f^) « 2 meV. The estimation for the holes tunneling parameter appears to 

be far less, taking rrihh = 0.5 mo, to'^,j = 0.15 toq from ([12]) one gets (t'^) ~ 0.01 meV. The polarization degree is to 
be estimated using ([^^ . We take A = 1 mcV, = Th = 20 K, the sheet concentration of the impurities n = 10^^ 
cm~^. Then for the case of the donor impurity t = , Eq ^ 4 meV, fih = —1 meV, /ie = , one gets \P\ « 40%, 
for the acceptor impurity t = , yUe = —1 meV, = 2 meV, = gives \P\ ~ 0.5%. An illustration of the 
luminescence spectra for the two circular polarizations is presented FigO For this we used an intermediate value for 
the tunnehng parameter = 0.3 meV (|P| w 0.15%) and accounted for inhomogeneous broadening of the spectra 
by normal distribution of the bandgap Eg with the dispersion cr = 3 meV (corresponds to the fluctuation of the QW 
width by half a monolayer). 

VII. SUMMARY 

The presented theory describes the tunnel coupling between a continuum of states in the QW and an impurity bound 
state located outside of the QW. We utilized the well known Fano approach for calculation of the matrix elements 
for the direct interband optical transitions in the QW. For such transitions the tunnel coupling of the 2D QW states 
with the impurity states leads to the drop of the luminescence spectral density at the frequency corresponding to 
the configuration resonance. This modification of the spectra leads to an integral circular polarization of the light 
emitted from the QW provided the bound hole state is split in the projection of the hole angular momentum. The key 
advantage of the approach used in the present study is that the unknown eigenfunctions of the system are expressed 
through those of the uncoupled states. Given the expansion any effects on the localized state can be translated 
into effects for the whole coupled system. For this reason it is capable of describing other effects expected in such 
systems like anisotropy of the holes g-factor in the QW induced by the paramagnetic impurity or the indirect exchange 
interaction between the bound states provided by the 2D free carriers in the QW. 
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8- layer 
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FIG. 1: (Color online) Mechanism of polarization of the luminescence for the acceptor type impurity. The localized hole levels 
split in magnetic field. Each of them effectively couples with the 2D holes having certain projection of angular momentum. 
Shifted positions of the resonances with account for temperature distribution of the holes cause the difference in intensities of 
circular polarizations cr^, a~ . The scheme also shows the simple electrostatic model described in the text. 
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FIG. 2: (Color online) Modification of the luminescence spectrum by tunneling configuration interaction. The integral polariza- 
tion occurs when the carriers distribution function (dashed line) strongly varies in the vicinity of the configuration resonances. 
cJo is the position of the resonance without bound level splitting. 
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FIG. 3: (Color online) An example of calculated luminescence spectra for the two circular polarizations. The case of anti- 
ferromagnetic coupling implies I~ corresponds to {o'^) polarization while /+ to (cr^) polarization. The parameters used in 
calculations are given in the text. 



